Multiband effects on /3-FeSe single crystals 
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We present the upper critical fields fj,oHc2{T) and Hall effect in /3-FeSe single crystals. The 
HoHc2{T) increases as the temperature is lowered for field applied parallel and perpendicular to (101), 
the natural growth facet of the crystal. The noHc2(T) for both field directions and the anisotropy 
at low temperature increase under pressure. Hole carriers are dominant at high magnetic fields. 
However, the contribution of electron-type carriers is significant at low fields and low temperature. 
Our results show that multiband effects dominate jj,oH C 2(T) and electronic transport in the normal 
state. 

PACS numbers: 74.70.Xa, 74.25.Op, 74.62.Fj 



I. INTRODUCTION 

The discovery of iron-based superconductors has gen- 
erated a great deal of interests because of rather high 
transition temperature T c and high upper critical fields 
11qH C 2- Crystal structures of iron-based superconduc- 
tors can be mainly categorized into several types: FePn- 
1111 type (REOFePn, RE = rare earth; Pn = P or As), 
FePn-122 type (AFe2As2, A = alkaline or alkaline-earth 
metals), FePnlll type (AFeAs), FeCh-11 type (FeCh, 
Ch = S, Se, Te), FeCh-122 type {A x Fe 2 -yCh 2 ), and 
other structures with more complex oxide layers De- 
spite structural similarity, i.e., shared FePn or FeCh 
tetrahedron layers, iron-based superconductors exhibit 
diverse physical properties. These include possible 
differences in pairing symmetry^ - — relation to com- 
peting or coexisting order states (spin density wave 
vs. superconductivity)^ - — and diverse normal state 
properties*^ FeCh-11 type materials are of special inter- 
est because their crystal structure has no blocking lay- 
ers in between FeCh layers, yet they have similar cal- 
culated Fermi surface topology when compared to other 
iron-based superconductors^ Furthermore, they also ex- 
hibit some exotic features, such as significant pressure 
effect ) 15 : 16 and excess Fe with local moment according 
to theoretical calculation^ 

The [1qH C 2 gives some important information on fun- 
damental superconducting properties: coherence length, 
anisotropy, details of underlying electronic structures 
and dimensionality of superconductivity as well as in- 
sights into the pair-breaking mechanism. Previous stud- 
ies on FeTei-aSez and FeTei-^Sa, single crystals indi- 
cate that the spin-paramagnetic effect is the main pair- 
breaking mechanism*^ - — However, for FePn-1111 and 
FePn-122 type superconductors the two-band effect with 
high diffusivity ratio between different bands dominates 
t, H c2 (T)2.-ZL 

On the other hand, magnetic penetration depth study 
of /3-FeSe polycrystal indicates that /3-FeSe is a two-band 
superconductor. 2 '* Therefore, it is of interest to inves- 



tigate multiband and spin paramagnetic effects on the 
lioH C 2 of /3-FeSe. An extremely complex binary alloy 
phase diagram and associated difficulties in single crys- 
tal preparation impeded the growth of pure /3-FeSe sin- 
gle crystals.— Hence, systematic studies of anisotropy in 
fiQH C 2(T) and pair breaking mechanism in high magnetic 
field are still lacking. 

In this work, we report on the upper critical fields of 
pure /3-FeSe single crystals in dc high magnetic fields up 
to 35 T at ambient and high pressures. The results shows 
that two-band features dominate the pair breaking with 
additional influence of spin paramagnetic effect. 



II. EXPERIMENT 

Details of crystal synthesis and characterization arc 
explained elsewhere^ 5 - The ^lqH C 2 is determined by mea- 
suring the magnetic field dependence of radio frequency 
(rf) contactless penetration depth in a static magnet up 
to 35 T at the National High Magnetic Field Labora- 
tory (NHMFL) in Tallahassee, Florida. The rf tech- 
nique is very sensitive to small changes in the rf pen- 
etration depth (about 1-5 nm) in the mixed state and 
thus is an accurate method for determining the [j.qH C 2 of 
superconductors fS 6 - At certain magnetic field, the probe 
detects the transition to the normal state by tracking the 
shift in resonant frequency, which is proportional to the 
change in penetration depth as A A oc AF. Small sin- 
gle crystals were chosen and the sample was placed in 
a circular detection coil. More details can be found in 
Refs. 29 and 30. For measurement under pressure, the 
sample was placed in a 15 turn coil within the gasket 
hole of a turnbuckle diamond anvil cell (DAC) made of 
beryllium copper and containing diamonds with 1.2 mm 
culets. 27 The pressure was calibrated at ^ 4 K by com- 
paring the fluorescence of a small chip of ruby within the 
DAC with an ambient ruby at the same temperature. 28 
The small dimensions of the DAC allow for angular rota- 
tion with respect to the applied magnetic field so H|| (101) 
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III. RESULTS AND DISCUSSION 
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FIG. 1. Field dependence of frequency shift (AF) for (a) 
H|| (101) and (b) H±(101) at various temperatures. AF(T 
— 10 K) is set as a normal-state background signal. Inset 
(a) and (b): enlarged parts near the points that deviate from 
background signals. poHc2(T) are determined from the in- 
terceptions between extrapolated straight lines below inflex- 
ion points and the curves of normal-state backgrounds for 
H 1 1 ( 101) and H_L(101), respectively. In the inset (a), the in- 
tersection of two slopes in AF(T — 6 K) curve exhibits an- 
other criterion to determine the poH C 2(T) and the difference 
between two criterions is taken as the error bar of poH C 2{T). 



and H_L(101) orientations can be explored in situ. Using 
four-probe configuration of Hall measurement, the Hall 
resistivity was extracted from the difference of transverse 
resistance measured at the positive and negative fields, 
i.e., p xy (H) = [p(+H) — p(—H)]/2, which can effectively 
eliminate the longitudinal resistivity component due to 
voltage probe misalignment. 



As shown in the main panel of Fig. 1(a) and (b), the 
rf shift (AF) at 10 K (above T c ) shows a smooth and al- 
most linear magnetic-field dependence without any steep 
changes for both field directions. In the normal state the 
rf shift is sensitive to the magnetoresistance of the sam- 
ple and detection coili^I However, when the temperature 
is below T c , there is a sudden increase of AF(H) which 
deviates from the background signal. This corresponds 
to entry to the mixed state. Moreover, with decreasing 
temperature, the inflexion points of the AF(H) curves 
shift to higher field for both field directions, consistent 
with the higher pqH C 2(T) at lower temperature. The 
temperature dependence of iiqH c i(T) for H|| (101) and 
H_L(101) is determined from the intersections of AF(H) 
curves between the extrapolated slopes of the rf signals 
below inflexion points and the normal-state backgrounds 
(T = 10 K) (insets in Fig. 1(a) and (b))M> The difference 
between this and other criterion, e.g. the intersection of 
extrapolated slopes below and above inflexion points in 
each AF(H) curve, is taken as the error bar of poH C 2{T). 

In order to compare the upper critical fields determined 
from different measurement methods, the poH C 2(T) ob- 
tained from AF(T,H) curves and p(T,H) data with 
different criteria are plotted together (Fig. 2(a)) M In 
the low field region, the temperature dependence of 
[loH&iT) determined from the rf shift is almost linear 
with slight upturn near T C (H — T). This is close to the 
HoH C 2,zero{T) determined from 10% p n (T,H). It is con- 
sistent with the results reported in the literature^ As- 
suming pqH C 2(T = 0.35 K) poH C 2(0), the zero temper- 
ature limit of upper critical fields are 17.4(2) and 19.7(4) 
T for H||(101) and H_L(101), respectively. On the other 
hand, according to the Werthamer-Helfand-Hohenberg 
(WHH) theory, orbital pair breaking field //oif C 2(0) = 
-0.693T c (dp H C 2/dT\T c )r^ an d using the initial slopes 
dpoH C 2/dT\T c at low fields obtained from p(T, H) data (- 
2.54(4) T/K for H||(101) and -2.55(4) T/K for H_L(101)) 
with T c = 8.7 we obtain the p H c2 (0) are 15.3(2) 
and 15.4(2) T for H||(101) and H_L(101), respectively. 
This is smaller than experimental results. The deviation 
from WHH model is clearly seen in Fig. 2(b), where the 
HoHcziT) becomes gradually larger than expected values 
from theory. The enhancement of the p>qH C 2 in the low 
temperature and high field region implies that multiband 
effect are not negligible. On the other hand, assuming the 
electron-phonon coupling parameter X e - P h = 0.5 (typi- 
cal value for weak-coupling BCS superconductors),— the 
Pauli limiting field (J. H P (0) = 1.86T c (l+A e _ p/l ) 1/2 is 19.8 
T<2£ This is nearly the same as the jUo-ff C 2,J?±(ioi) (T — 
0.35 K) and larger than values for H|| (101) or the or- 
bital pair breaking fields. It suggests that the spin- 
paramagnetic effect might also have some influence on 
the upper critical fields. This is rather different from 
other FeCh-11 superconductors where the Pauli limit- 
ing fields are much smaller than orbital pair-breaking 
fields and therefore the spin-paramagnetic effect governs 



3 




0123456789 10 



T(K) 

FIG. 2. (a) Temperature dependence of poH C 2{T) of /3- 
FeSe single crystal for H|| (101) (closed symbols) and H_L(101) 
(open symbols) obtained from p(T) and AF curves, (b) Fits 
of HoH C 2(T) for H|| (101) using eq. (1) for different pairing 
scenarios: (1) WHH; (2) vj > 0, An = 0.241, A 22 = 0.195, 
A12 = A21 = 0.01, rj = D2/D1 = 0.40; (3) w > 0, An = A 22 
= 0.5, A12 = A21 = 0.25, r, = D2/D1 = 0.44; (4) w > 0, An 
= 0.8 A22 = 0.34, A12 = A21 = 0.18, rj = D2/D1 = 0.32; and 
(5) vj < 0, An = A22 = 0.49, A12 = A21 = 0.5, rj = D2/D1 = 
0.35; 

According to the two-band BCS model in the dirty 
limit with orbital pair breaking and negligible interband 
scattering, 3 — [i^Hci is given by 

a [lni + U(h)][\nt + U(r)h)] + a 2 [lrrt + U(r)h)] 

+ a 1 [lnt + U(h)]=0 (1) 

where t = T/T c , U(x) = tp(l/2 + x) - ijj(x), ip(x) the 
digamma function, i] = D^j ' D\, D\ and D 2 are intraband 
diffusivities of the bands 1 and 2, h = H c2 D\j (2<j) T), 
4>q the magnetic flux quantum, ao, a%, and a 2 are 
constants described with intraband- and interband cou- 
pling strength, ao = 2m /Ao, a% = 1 + A_/Ao, and 
oi = 1 — A_/A , where w = A11A22 - Ai 2 A 2 i, A = 
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FIG. 3. (a) Field dependence of p xy (H) at various temper- 
atures. Solid lines are the fitting results using eq. (2). for 
T < 60 K and single-band model for T = 60 K. In order to 
exhibit data clearly, the p xy (H) at different temperatures are 
shifted along vertical axis with certain values, (b) Temper- 
ature dependence of carrier density n(= nn — n e ) of /3-FeSe 
crystal. 



(A^ + 4.A12A21) 1 / 2 , and A_ = An - A 2 2- Terms An and 
A22 are the intraband couplings in the bands 1 and 2 
and A12 and A21 describe the interband couplings be- 
tween bands 1 and 2. It should be noted that if rj = 
1, eq. (1) will reduce to the simplified WHH equa- 
tion for single-band dirty superconductors^ By using 
the coupling constants determined from an riSR exper- 
iment with very small interband coupling^ the com- 
bined P-qH c2 ^h\\{ioi){T) data from both rf and resistiv- 
ity measurements can be very well explained (Fig. 2(a) 
fit lines). The ratio of band diffusivities is 77 = 0.40, 
which is similar to the value of FeAs-122 but much larger 
than that of other two-band iron-based superconduc- 
tors, such as FeAs-lllli2i - — With current coupling con- 
stants, it leads to the similar shape of A t o-ff C 2,ff||(ioi) (T) 
when compared to the FeAs-122^ but significantly dif- 
ferent from FeAs-1111 where there is an obvious upturn 
at low temperature<2i^2 We have also performed fits for 
different values of coupling constants^ - — (1) dominant 
intraband coupling, vj > and (2) dominant interband 
coupling, w < 0. The different sets of fitting parameters 
result in almost identical result, fitting the experimental 
data well (Fig. 2(b)). The derived 77 is in the range of 
0.32-0.44, suggesting that the fitting results are insen- 
sitive to the choice of coupling constants. Thus, either 
interband and intraband coupling strength are compa- 
rable or their difference is below the resolution of our 
experiment. 

In order to further investigate multiband character- 
istics in /3-FeSe, we studied the Hall effect of /3-FeSe 
(Fig. 3). According to the band calculations, at least 
four bands originated from Fe 3d orbitals cross the Fermi 
levelJ^ Two bands are hole type and the other two are 
electron typeJ^ We use a simplified two-carrier model 
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including one electron type with electron density n e and 
mobility p e and one hole type with hole density n-h and 
mobility ph- According to the classical expression for 
the Hall coefficient including both electron and hole type 
carriers^ 



R _ 1 {pjn h - p 2 n e ) + (p h p e ) 2 (p H) 2 (n h - n e ) 
e (p e n h + p h n e ) 2 + (p h p e ) 2 (p H) 2 (n h - n e ) 2 

(2) 

Once there are two carrier types present, the field de- 
pendence of p xy {H) will become nonlinear. Moreover, 
eq. (2) gives R H = e _1 • (p 2 h n h - p? e n e )/ (Me^fc + Ph,n e ) 2 
when hqH — > 0, and Rh = e _1 -l/(n/ l — n e ) when pgH — >■ 
oo. As shown in inset (a) of Fig. 3, p xy {H) is positive 
and almost linear in poH at T = 60 K, indicating the 
hole type carrier is dominant. However, p xy (H) exhibits 
obvious nonlinear behavior below 50 K and even changes 
sign in low fields at 15 K (inset (b) in Fig. 3). This is a 
signature of coexistence of electron and hole type carri- 
ers. The p X y(H) can be described very well using a linear 
relation for T — 60 K and eq. (2) for T ^ 50 K as shown 
with the solid fit lines in the inset (a) and (b) of Fig. 

3. The obtained carrier density n(= Uh — n e ) changes 
from 1.93xl0 21 cm" 3 (15 K) to 4.7xl0 21 cm" 3 (60 K) 
gradually. The change of sign of p xy (H) in the low field 
region at 15 K indicates {p^rih — p^n e ) < 0. Because 
nh — n e > at higher field, it indicates that the p e > pt 
at low temperature, consistent with the band structure 
calculation results.— Moreover, the negative Seebeck co- 
efficients in /3-FeSe below ~ 250 K also confirm that the 
electron band is dominant at low temperature^ 

Since there is remarkable pressure effect on T c for /3- 
FeSer^J^ it is instructive to study the pressure depen- 
dence of upper critical fields. As shown in the inset of Fig. 

4, under pressure (P — 0.51 GPa), the inflexion point of 
AF(H) curve shifts to higher field when compared to the 
ambient pressure curve, suggesting that the pjQH&iT) is 
enhanced with pressure. It is consistent with the signif- 
icant positive pressure effect of T c for /3-FeSeJ^& The 
temperature dependence of PqH c2 (T) for H||(101) and 
H_L(101) shows that the upper critical fields for both field 
directions are enhanced in the whole measured tempera- 
ture region under pressure. The poH C 2(T = 0.45 K) for 
H_L(101) is about 24.6 T, close to the estimated value 
at 1.48 GPa using linear extrapolation^ It suggests 
that poH c2 (Q) at 0.51 GPa should be larger than linear- 
extrapolated value. This could originate from the dif- 
ference in sample purity between our single crystals and 
polycrystals or intrinsic multiband effect. As shown in the 
inset (b) of Fig. 4, at ambient pressure, the anisotropy 
of p a H c2 (T), 7 (P = GPa, T) = p Q H c2>H±{wl) {P = 
GP&,T)/p H c2 , Hmol) (P = GPa, T), is smaller than 
in other iron based superconductors, especially at high 
temperature. Moreover, the temperature dependence 
of j(P, T) increases at high temperature and decreases 
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FIG. 4. (a) Temperature dependence of p H c2 (T) for H||(101) 
(closed symbols) and H_L(101) (open symbols) at ambient 
pressure and P — 0.51 GPa obtained from AF curves. In- 
set (a) field dependence of AF at and 0.51 GPa for H||(101) 
at T = 0.35 K. Inset (b) The temperature dependence of the 
anisotropy of poHc2{P,T) at P = and 0.51 Gpa. 



when T <C T c , which is different from other iron based 
superconductors in which the j(T) usually decreases with 
temperature. The increase of j(P, T) with tempera- 
ture has also been observed in two-band superconduc- 
tor MgB2. This may be due to the higher contribution 
of the band with lower band anisotropy at low temper- 
ature. The decrease of j(P,T) with temperature when 
temperature is far from Tc may be related to the possi- 
ble spin-paramagnetic effect. On the other hand, under 
pressure, the j(P — 0.51 GPa, T) increases when com- 
pared to the value at ambient pressure. This could orig- 
inate from the pressure-induced Fermi surface changes 
that increase the anisotropy of Fermi velocity (diffusiv- 
ity) of dominant band. It should be noted that in order to 
study the anisotropy and pressure evolution of pqH C 2(T) 
more clearly, the pressure dependence of PqH c2 (T) along 
crystallographic axes should be measured in the future. 



IV. CONCLUSION 

In summary, we studied the upper critical field of f3- 
FeSe crystals. The results indicate that the two band 
effects dominate the poH c2 (T), with possible influence of 
spin-paramagnetic effect. A nonlinear field dependence 
of p xy (H) at low temperature also confirms the existence 
of multiple bands in electronic transport. The dominant 
carriers are hole-type in high field but electron type carri- 
ers become important in low field due to either increased 
carrier density or enhanced mobility. The PqH c2 (T) is 
enhanced for both field directions and the anisotropy of 
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IM)H C 2{fy is also increased under pressure. 
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